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Given a stable Levy process X = {Xt)a<t<T of index a G (1,2) 
with no negative jumps, and letting St = supp^^^j Xs denote its run- 
ning supremum for t G [0, T] , we consider the optimal prediction prob- 
lem 

inf E(ST-X^f, 

0<T<T 

where the infimum is taken over all stopping times r of X, and 
the error parameter p £ (1,0;) is given and fixed. Reducing the op- 
timal prediction problem to a fractional free-boundary problem of 
Riemann-Liouville type, and finding an explicit solution to the lat- 
ter, we show that there exists a, G (1,2) (equal to 1.57 approxi- 
mately) and a strictly increasing function p*:(a,,2) — >■ (1,2) satis- 
fying p,(a*-f) = 1, p,(2— ) — 2 and p*(a) < a for a £ (a*, 2) such 
that for every a £ (a«,2) and p G (l,p,(a)) the following stopping 
time is optimal 

T,=inf{te[0,T]:St-Xt>z4T-t)'-^°'}, 

where 2:* £ (0, 00) is the unique root to a transcendental equation 
(with parameters a and p). Moreover, if either aG(l,a*) orpG 
{p,{a),a) then it is not optimal to stop at t G [0,r) when St — Xt is 
sufficiently large. The existence of the breakdown points a, and p. (a) 
stands in sharp contrast with the Brownian motion case (formally cor- 
responding to a — 2), and the phenomenon itself may be attributed 
to the interplay between the jump structure (admitting a transition 
from lighter to heavier tails) and the individual preferences (repre- 
sented by the error parameter p). 
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1. Introduction. Stopping a stocliastic process X = {Xt)o<t<T as close 
as possible to its ultimate supremum St = sviY>Q<s<T^s is an objective of 
both practical and theoretical interest. Speaking in general terms, the opti- 
mal prediction problem can be formulated as follows 

(1.1) V= inf diXr^Sr), 

0<T<T 

where the infimum is taken over all stopping times r of X, and d is a dis- 
tance/error function [e.g., d{Xr,ST) = E(5t — Xr)^ where p > is a param- 
eter quantifying the error]. Variants of these problems have been studied in 
the past mostly in discrete time (see, e.g., [4, 9, 12, 14]), and the case of 
continuous time has been studied in the recent papers [11] and [17] when X 
is a standard Brownian motion. This study was extended in [6] to the case 
of Brownian motion with drift. It was observed there that the existence of 
a nonzero drift leads to optimal stopping boundaries having a complex struc- 
ture which in some cases appears to be counter-intuitive. For other optimal 
prediction problems studied to date, we refer to [7, 8, 22-24] (see also [18], 
Chapter VIII). In these problems, it is assumed that the underlying process 
has continuous sample paths. 

The purpose of the present paper is to initiate a study of the optimal 
prediction problems for processes with jumps in continuous time, and to 
examine the extent to which the jump structure influences the resulting op- 
timal stopping boundaries. To stay close to the more familiar case of Brow- 
nian motion, we study the case when X is a stable Levy process of index 
a € (1,2), and to focus on one particular aspect of the jump structure we 
consider the case when X jumps upward only (i.e., when X has no negative 
jumps). It turns out that already these hypotheses lead to a complicated 
optimal prediction problem, which apart from initial similarities with the 
case of Brownian motion (through the scaling property and deterministic 
time-change arguments) requires novel arguments to be developed in order 
to find a solution. These complications are primarily attributed to the un- 
derlying jump structure which leads to the relatively unexplored avenue of 
integro-differential equations (fractional calculus) instead of more familiar 
differential equations. Yet another difficulty (that the law of St was not 
available in the literature prior to the present study) is now overcome by 
the accompanying paper [2], and the knowledge of this law plays a key role 
in our treatment of the optimal prediction problem below. 

Our main findings (Theorem 11) can be summarized as follows. Given 
a stable Levy process X = {Xt)o<t<T of index a G (1,2) with no negative 
jumps, and letting St = supQ^^^^Xg denote its running supremum for t G 
[0,T], we consider the optimal prediction problem 



(1.2) 



V= inf EiST-XrY, 

0<T<T 
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where the infimum is taken over all stopping times r of X, and the er- 
ror parameter p E (1,0;) is given and fixed (we will see in Section 2 below 
why the restriction to this interval is natural). Reducing the optimal pre- 
diction problem to a fractional free-boundary problem of Riemann-Liouville 
type, and finding an explicit solution to the latter, we show that there ex- 
ists a* € (1,2) (equal to 1.57 approximately) and a strictly increasing func- 
tion : (a*, 2) — )• (1, 2) satisfying p*(q;*+) = l,p*(2-) = 2 and p*{a) < a for 
a G (a*, 2) such that for every a € (a=i,,2) and p £ (l,p*(a)) the following 
stopping time is optimal 



where G (0, oo) is the unique root to a transcendental equation (with pa- 
rameters a and p) . This extends the analogous results for a standard Brown- 
ian motion X derived in [11] and [17] when p = 2 and p G (1,2), respectively. 
Moreover, if either a G (1, a*) or p G (p*(a), a) then it is not optimal to stop 
at i G [0, T) when St — Xt is sufficiently large. The existence of the breakdown 
points a=f and (a) stands in sharp contrast with the Brownian motion case 
(formally corresponding to a = 2), and the phenomenon itself may be at- 
tributed to the interplay between the jump structure (admitting a transition 
from lighter to heavier tails) and the individual preferences (represented by 
the error parameter p). In particular, recalling that the index a quantifies 
the heaviness of the upward tails of the process X, we see that the result 
may be broadly interpreted as follows: the heavier the upward tails the larger 
the optimal stopping time. While this conclusion is close to naive intuition, 
and the interpretation itself may also be extended to account for the individ- 
ual preferences, the fact that the solution method can detect the breakdown 
points exactly appears to be of considerable practical and theoretical inter- 
est. Other interesting features of the problem include the remarkable proba- 
bilistic representation of the solution to the Ito/Riemann-Liouville/Caputo 
free-boundary problem that is novel in the case of Brownian motion as well. 

2. The optimal prediction problem. 1. Let X = {Xt)t>o be a stable Levy 
process of index a G (1,2) whose characteristic function is given by 



for A G M and t>0 with c > 0. Let S = {St)t>o denote the supremum process 
of X, that is. 



(1.3) 



n = inf{t G [0, T]:St-Xt> z^{T - t)^/"} 



(2.1) 




(2.2) St = sup Xs 



0<s<t 



for t > 0. Consider the optimal prediction problem 



2.3) V= inf EiSr-XrY, 



Q<T<T 
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where the infimum is taken over all stopping times r of X [i.e., stopping 
times with respect to the natural filtration = (t{Xs : < s < t) generated 
by X for i > 0] . It is assumed in (2.3) that the error parameter p € (1, a) and 
the terminal time T > are given and fixed (we will see below that there is 
no restriction in assuming that T = 1). 

2. The following properties of X are readily deduced from (2.1) using 
standard means (see, e.g., [3] and [15]): the law of {X^t)t>o is the same as 
the law of (a^/^Xt) t>o for each cr > given and fixed (scaling property); 
X is a martingale with EXt = for all t > 0; X jumps upward (only) and 
creeps downward [in the sense that P{Xp^ = x) = 1 for x < where px = 
inf{t > 0:^4 < x} is the first entry time of X into (— oo,x)]; X has sample 
paths of unbounded variation; X oscillates from — oo to +oo (in the sense 
that lim inf t_).(X) = — oo and lim sup^^o^ = +oo both a.s.); the starting 
point of X is regular [for both (— oo,0) and (0,+oo)]. Note also that the 
Levy measure v oi X equals 

(2.4) T,{^dx) = ^dx 

on the Borel u-algebra of (0,oo). Setting, for example, c= l/(2r(— a)) we 
see from (2.1) that X = X[a) converges in law to a standard Brownian 
motion i? as at 2- We moreover see from (2.4) that when a is closer to 2 
then the (upward) jumps of X have lighter tails, and when a is closer to 1 
then the (upward) jumps of X have heavier tails. Thus, in many ways, the 
process X resembles a standard Brownian motion B, however, the existence 
of (upward) jumps of X represents a notable exception. Note also that Xt 
is not equal in law to —Xt for fixed t > unlike in the case of B. 

3. The error parameter p in the problem (2.3) is assumed to belong to 
(l,a) for two reasons. First, it is well known (see, e.g., [21], page 159) 
that for a Levy process X = {Xt)t>{) and a number p > given and fixed, 
the following three facts are equivalent: (i) EXf < oo for some/all i > 0; 

(ii) Esupo<5<(Xf < oo for some/all t > 0; (iii) /-^ x^v{dx) < oo. In the case 
of our process X when is given by (2.4) above, it is easily seen that 

(iii) holds [and thus both expected values in (i) and (ii) are finite] if and 
only if p < a. In particular, the latter condition then also implies that the 
value V in (2.3) is finite. Second, if p = 1 then the optimal prediction prob- 
lem (2.3) is trivial since EXr = for every (bounded) stopping time r of X 
due to the martingale property of X. Hence, p € (l,a) represents a natural 
assumption on the error parameter. 

4. Note that there is no loss of generality if we assume that T = 1 in the 
problem (2.3). Indeed, if we set V = V{T) to indicate dependence on T > 
in (2.3), then by the scaling property of X we see that V{T) = Tp/°'V{1) and 
there is a simple one-to-one correspondence between the stopping times r in 
the problem V{T) and the stopping times a in the problem V{1) (obtained 
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by setting a = t/T). For this reason, we wih often assume in the sequel that 
the horizon T in (2.3) equals 1. 

5. Projecting future onto present. One of the key initial difficulties in the 
optimal prediction problem (2.3) is that the expression after the expectation 
sign contains the random variable St and as such depends on the (ultimate) 
future of the process X that is unknown at the present (stopping) time r € 
[0,T). In our first step therefore (similarly to [11] and [17]), we will project 
the future states of X onto the present /past states of X by conditioning with 
respect to J-^ and exploiting stationary/independent increments oi X. As 
already mentioned above, we may and do assume that T = 1 in the sequel. 

To this end, note that we have 



E{{Si-xtm 



(2.5) 



;(( sup {Xs - Xt) V sup {Xs - X^)) "ijf 

^ ^0<s<t t<s<l ' 



(E(2/V5: 



i-i 



\y=St-Xt 



since supi<^<i(X^ - Xt 



S\-t is independent from and St 



Xt is 



/"/''-measurable. Moreover, we can write 



E(yV5: 



(2.6) 



/>oo 

: / ?({yy S^.tY>z)dz 
Jo 

/•oo 

P{SU>z)dz 

JyP 

POO 

yP+ P{{l-t)P/''Sf>z)dz 

JyP 

(i-*)"°[((t4w)' 

fOO 

P{SP>w)dw 



(y/(l-t)V-)P 



--■Fit,y) 



upon using that Si-t ^= (1 — tY^^'Si by the scaling property of X and sub- 
stituting w = z/{l - t)P/°'. Combining (2.5) and (2.6), we get 



(2.7) 



E((5i 



XtTlTt' 



F{t,St-Xt) 



for all t > 0. Using the fact that each stopping time r of X is the limit 
of a decreasing sequence of discrete stopping times r„ of X as n — )• oo, it 
is easily verified using Hunt's lemma (see, e.g., [26], page 236) that (2.7) 
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extends as follows 

(2.8) E{{Sl-Xrr\T^)=F{T,Sr-Xr) 

for all stopping times t oi X with values in [0, 1]. Setting 

(2.9) Yt = St- Xt 

for t > it is well known (see, e.g., [3]) that Y = (Yt)t>o is a time-homogeneous 
(strong) Markov process with respect to {J^i^)t>o (obtained by reflecting X 
at its supremum S). Taking E on both sides in (2.8) and using the nota- 
tion (2.9), we see that the optimal prediction problem (2.3) reduces to the 
optimal stopping problem 

(2.10) V= inf EF(T,Yr), 

0<T<1 

where the infimum is taken over all stopping times t X. This optimal 
stopping problem is two-dimensional (see, e.g., [18], Section 6) since the 
underlying (strong) Markov process is the time-space process ((t, lt))o<t<i 
and the horizon 1 is finite. We will now show (similarly to [11]) that this 
problem can further be reduced to a one-dimensional infinite-horizon optimal 
stopping problem for a (killed) Markov process Z = (Zs)s>o- It should be 
noted that the time-change arguments used in [11] when X is a standard 
Brownian motion are not directly applicable in the present context (due to 
the absence of Levy's characterization theorem). 

5. Deterministic time change. Motivated by the form of the function F 
in (2.6), we now introduce the deterministic time change 

(2.11) ^(s) = l-e-"^ 

where t{s) € [0, 1) is the "old" time and s € [0, oo) is a "new" time. Note 
that r = t{a) is a stopping time with respect to {T^)t>o if and only if 
a = t^~^\T) is a stopping time with respect to {J^^g-^)s>o- Letting denote 
the distribution function of Si and setting 

(2.12) Giz) = E{zV Si)P = zP+ / {1 - Fs.iw^^P)) dw 

for z > 0, we see from (2.6) and (2.12) that 

(2.13) F{t,St-Xt) = e-P'G{Zs) 

for all t = t{s) € [0, 1) and all s € [0, oo) satisfying (2.11), where Z = (Zs)s>o 
is a new stochastic process defined by 

(2.14) Z, = e^(%)-Xi(,)) 

for s > 0. It turns out that Z is a time- homogeneous (strong) Markov pro- 
cess. Moreover, the following proposition reveals that one can enable Z to 
start at arbitrary points and still preserve the (strong) Markov property. 
This fact will play a prominent role in the main proof below. 
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Proposition 1. The stochastic process Z = {Zs)s>o defined in (2.14) 
is a time-homogenous (strong) Markov process with respect to the filtration 
{J^^g))s>0- Moreover, if we set 

(2.15) = e^(zV%)-Xi(,)) 

for s>Q and z S M+, then Pz ■= Law((Z|)s>o|P) defines a family of prob- 
ability measures on the canonical space of cadldg functions {D^,B{D^)) 
under which the coordinate process C = {Cs)s>o is (strong) Markov with 
Pz{Co = z) = l for ze^+. 

Proof. We have 



^s+h 



^s+h 



{z V S't(j,+/i) - 



,s+h 



{Z V St{s 



(2.16) 



■^t(s))V( sup {Xr-Xt(^sy 
H(s)<r<t(s+h) 



Zf V e" I sup 

t(s)<r<t(s+h) 



{Xr - Xt(^s))^ - e''{Xt(s_^,h) - ^t{s))^ 



for s > and h>0 given and fixed. By stationary independent increments 
and the scaling property of X, we see that 



sup (Xr-Xt(^)) 

t{s)<r<t{s+h) 



sup {Xr — Xi_^-a3) 



law 



(2.17) 



sup Xr 

0<r<e-"»(l-e-«'») 



law 



sup v\(-^gas-)/g. 

0<re°"=<l-e-"'' 



law 



e sup yi-r 

0<r-<l-e-"'i 



Xr = e-'S, 



and likewise 
(2.18) 



^t{s+h) — ^t{s) 



^1- 



,ga(s + /i) 



X, 



law 



(l_e-ah) 



law 



both being independent from ^^^^.y Combining (2.16)-(2.18), we get 
(2.19) E(/(Z,^+,)|J-^,)) = E(/(e'^(u;V5i(,)-X,(,))))U^. 

for any (bounded) measurable function / : IR_|_ M from where ah the claims 
follow by standard means [observe that the deterministic function on the 
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right-hand side of (2.19) does not depend on s (implying that Z is a time- 
homogenous Markov process) as well as that it defines a continuous and 
bounded function of w whenever / is so (Feller property) implying that Z 
is a strong Markov process]. This completes the proof. □ 

Note from (2.14) that Z is a transient process (satisfying Zs ^ oo as 
s — 7> oo) having downward jumps only (since X jumps upward). The state 
space of Z equals 1R+. 

3. The optimal stopping problem. 1. From (2.10) and (2.13), we see 
that the optimal prediction problem (2.3) reduces to the optimal stopping 
problem 

(3.1) V= inf Ee-P^GiZ^), 

0<o-<oo 

where the infimum is taken over all stopping times a with respect to {^t^g^)s>o- 
This optimal stopping problem is one-dimensional and the horizon is infi- 
nite. The exponential term {e~P^)s>o in (3.1) corresponds to a new (strong) 
Markov process Z which may be identified with Z killed at rate p. 

2. To tackle the problem (3.1), we need to enable Z to start at any point 
in the state space IR+. This can be done using the result of Proposition 1 
above, and it leads to the following variational extension of (3.1): 

(3.2) V{z)= inf E,e-P'^G(Z^), 

0<(T<OO 

where the infimum is taken over all stopping times a with respect to {^t^s^)s>Oj 
and the process Z starts at z under P^. Moreover, by the result of Proposi- 
tion 1 we know that can be realized by (2.15) in terms of Z^ = {Zg)s>o 
under P, and this fact will be useful below when analysing properties of the 
mapping z i— )■ V{z) on M4.. 

3. Before we turn to a more detailed analysis of the problem (3.2), let 
us state some basic properties of G and V that will be useful throughout. 
Recall that f{z) ~ g{z) as z ^ zq means that liuiz^zQ f{z)/g{z) = 1 for zq G 
[—00, 00]. 

Proposition 2. The gain function G from (2.12) above and the value 
function V from (3.2) above satisfy the following properties: 

z^G{z) is (strictly) increasing and convex on 
with G(0) = E5[ > 0; 

z I— )• V{z) is increasing and continuous on M4.; 

(3.5) zP <V{z)<G{z) for all z€R+; 

(3.6) G{z) ~ zP and V{z) ^z^ as z^ 00. 



(3.3) 
(3.4) 
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Proof. Equation (3.3): recalling that Fs-^ denotes the distribution func- 
tion of Si, and letting fs^ denote the density function of 5*1, we find from 
the final expression in (2.12) that G'{z) = pz^-'^Fs^iz) > and G"{z) = 
p{p— 1) X 2:P~^F5^(2;) +pzP~^fs-i^{z) > for ah 2 > implying that z i-> G{z) 
is (strictly) increasing and convex, respectively. Likewise, we also see from 
the middle expression in (2.12) that G(0) = ES^ > as claimed. 

Equation (3.4): letting o" be a given and fixed stopping time, we see 
from (2.15) that z i— )■ Z| is increasing so that z i— )• G{Z^) is increasing, and 
the fact that z i— > V{z) is increasing follows directly from the definition (3.2). 
To show that z ^ y{z) is continuous, take zi < Z2 in M+ and note by the 
mean value theorem and (2.15) that 

< G(Z^^) - G(Z^O = G'mzi- - 

(3.7) =G'(0e'^(^2V5t(,)-ziV5t(,)) 

where i G (Zf ,Z|2). Since < ^ < e'^(z2 V 5i - /i), where we set h = 
inffxtKi Xf, it follows from (3.7) that 

(3.8) < Ee-P^GiZ^') - Ee^^'^ G{Zl^) <pE{z2 V Si - hY~\z2 - zi). 
Taking the infimum over all stopping times cr it follows that 

(3.9) ^<V{z2)-V{zi)<K{z2-zi), 

where K = pE{z2 V Si— IiY~^ < oo. This implies that V is continuous on 
(as well as Lipschitz continuous on compact sets in M_|_). 

Equation (3.5): the second inequality is obvious so let us derive the first 
inequality. For this, fix any z S M+ and note that G{z) > z^ and Jensen's 
inequality imply that 

l-(z)> inf Ee-^"^(Z^)^>( inf Ee~- Z^Y 

0<o-<oo \0<(T<OO / 

(3.10) =( inf E(zV5,(,)-X,(,)))'' 

\0<cr<OO / 

= ( inf E(zV Sr-Xr)Y = zP 

\0<r<l / 

upon using that there is a one-to-one correspondence between a and r as 
stated following (2.11) above. Note also that for the final equality we use 
the fact that EXr = since X is a martingale. This establishes the first 
inequality in (3.5) as claimed. 

Equation (3.6): note that (2.12) above implies that G{z)/ z'^ — )■ 1 as z — )■ oo, 
so that V{z)/ z'P — )• 1 as z — >■ oo follows by (3.5). This completes the proof. 
□ 
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4. Existence of an optimal stopping time. General theory of optimal stop- 
ping for Markov processes (see, e.g., [18]) can be used to establish the 
existence of an optimal stopping time in the problem (3.2). For this, let 
C = {z £ W+:V{z) < G{z)} denote the (open) continuation set, let D = 
{z e : V{z) = G{z)} denote the (closed) stopping set, and note that 

(3.11) E(supe-P'G{ZI)) <oo 

since e-P^G(Z|) = e-P^((Zf )p + P(5^ > w) dw) < (z V 5i - Ii)p + E5f 

for all s > 0, and the latter random variable clearly is integrable for each z € 
IR4.. Moreover, by (3.3) and (3.4) we know that the gain function z 1— )■ G{z) 
is lower semicontinuous on R_|_ and the value function z ^V{z) is upper 
semicontinuous on R_|_. Hence, by Corollary 2.9 and Remark 2.10 in [18], 
pages 46-48, we can conclude that the first entry time of Z into D given by 

(3.12) o-z) = inf {s > : G £)} 

is an optimal stopping time in (3.2). This stopping time is not necessarily 
finite valued [when the set in (3.12) is empty] and the value e~*"^°G(Z|^) 
in (3.2) can be formally assigned as {z\/ Si — XiY when gd = 00 since 
by (2.12) and (2.15) we have 

(3.13) e-P'G{Z^^)^{zy Si-XiY 

as s — )■ 00. This is in agreement with the usual hypothesis from general 
theory introduced to cover the case of infinite- valued stopping times. 

5. In addition to these general facts, it may be noted that the optimal 
stopping problem (3.2) plays an auxiliary role in tackling the optimal pre- 
diction problem (2.3), and it is clear from our considerations above that we 
only need to compute V{z) for z = 0. Thus, if we set = inf D then either 

< 00 when L> / (so that z^ £ D since D is closed) or z^, = 00 when 
D = 0. In the first case (when D 7^ 0), the first entry time of Z to 2* given 
by 

(3.14) 0-^, = inf{s >0:Zs = z^} 

is optimal in (3.2) under for z = 0. It should be recalled here that Z 
jumps downward only and creeps upward in so that Z will hit any point 
in (0, 00) with probability one due to its transience to +00. Recalling further 
the time change (2.11) we see that (3.14) translates into the fact that the 
stopping time 

(3.15) n = inf{t G [0, l]:St-Xt> z^{l - t)^/"} 

is optimal in (2.3) with T = 1. In the second case (when D = 0), we see 
that the optimal stopping time a^^. in (3.2) equals +00 under for z = 0. 
In this case, we have 

(3.16) Viz) = E{zV Si-Xif 
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for all z € M_|_ and the time change (2.11) implies that = 1 is optimal 
in (2.3) with 

(3.17) T^ = E(5i-Xi)P. 

A central question therefore becomes to examine when [0, 2^,) C C with 
z-f ^ D (it will be shown in Section 5 below that z^ cannot be zero). We 
will tackle this question by forming a free-boundary problem on [0, z^) for V 
defined in (3.2). For this, we first need to determine the infinitesimal char- 
acteristics of Z. 

4. The free-boundary problem. 1. The following proposition determines 
the action of the infinitesimal generator of the process Z defined in (2.14) 
in terms of the action of the infinitesimal generator of the reflected process 
Y = S — X . Below we let C|(M+) denote the class of twice continuously 
differentiable functions F : M_|_ — >■ R such that F' and F" are bounded on M+ . 

Proposition 3. The infinitesimal generator of the process Z is 
given by 

(4.1) l.zF{z) = zF'{z) + al.YF{z) 

for any F € C^(]R+) satisfying (4-6) below, where hy denotes the infinitesi- 
mal generator of the process Y . 

Proof. By the mean value theorem, we have 

lLzF{z)=lun-E{FiZl)-Fiz)) 

slO S 

= hm iE(F(e^(z V %) - X^^,^)) - F{z V - X,(,)) 

(4.2) +F(zV%)-X,(,))-F(z)) 
= hm ^E(F'(6)(z V Stis) - ^t(.))) 

= zF'(z) + aLyF(z), 

where for the second last limit we use that (e^ — l)/s — )■ 1 and F'{^s) F"{z) 
as s i since S (z V S'^j-g) — Xf-(^g^,e^{z V 5't(^) — X^j-^))), and for the last limit 
we use that t{s)/s — t- a as s J, and the result of Proposition 4 below. This 
completes the proof. □ 

2. The following proposition determines the action of the infinitesimal 
generator of the refiected process Y = S — X . We refer to the Appendix for 
the analogous result in the case of a general (strictly) stable Levy process X. 
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Proposition 4. The infinitesimal generator Ly of the reflected process 
Y = S — X takes any of the following three forms for y > given and fixed: 

ltd 's form 

ry 

(4.3) hYF{y)= {F{y-x)-F{y)+F'{y)x)^dx 

Jo X 

, c(F(0)-F(y)) ^ cF'{y) 



ay"' {a — l)y" ^ ' 

Riemann-Liouville 's form 

TT N c d^ fy Fix) ^ cF(0) 

(4.4) LyF(y) = — r-r^ / 7 \ ^ , dx ^ —, 

^ ' ^ ' a(a - 1) dy^ Jo {y - x)"-i ay" 

Caputo 's form 

c I'y F"(x) 

(4.5) LyF(y) = — ^ / dx, 

a{a-l) Jq (y-x)" ^ 

whenever F £ C^(M+) satisfies 

(4.6) F'(0+) = (normal reflection). 

Proof. It is enough to establish (4.3) since (4.4) and (4.5) can then 
be derived by (repeated) integration by parts using (4.6) (note that the 
equivalence of (4.3)-(4.5) under (4.6) remain valid for any F € C^[0, 00) PI 
C^(0, 00) satisfying \F"{x)\ = 0(x"~^) as x^O since a — 2 > -1). For this, 
fix t > and note that by Ito's formula we have 

F{Yt) = F{Yo)+ f F'{Ys-)dY, 

(4.7) 

+ ^ (F(n)-F(n„)-nn_)Ay,) 

0<s<t 

since [y, = 0. Indeed, the latter equality follows by recalling that X is 
a quadratic pure jump semimartingale (i.e., [X, X]'^ = 0) since it is a Levy 
process with no Brownian component (see [19], page 71), the process S is 
a quadratic pure jump semimartingale since it is of bounded variation (see 
Theorem 26 in [19], page 71), and the sum/difference of two quadratic pure 
jump semimartingales is a quadratic pure jump semimartingale (this can be 
easily verified using Theorem 28 in [19], page 75, e.g.). 

Since X jumps upward and creeps downward, it follows that dSs = ASg 
in terms of a suggestive notation, and hence from (4.7) we get 

(4.8) F{Yt) = F{Yo) + Mt+ ^ (F(y,_ + A^) - F(n_) + F'(y,_)AX,), 

0<s<t 
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where Mf = — J^* F'{Ys-) dXg is a local martingale for t > 0. By the BDG in- 
equality (see, e.g., [18], page 63) combined with the facts that F' is bounded 
on R+ and E[X, X]^ < oo with q = 1/2 since [X, X] is a stable process of 
index a/2> q [with Levy measure cdx / {2x^^°^/'^) as is easily verified di- 
rectly from definition] it follows that Esupo<s<j|Ms| < oo and hence M is 
a martingale. The right-hand side of this identity can be further rewritten 
as follows 

F{Yt) = F{Yo) + Mt 



+ > ([F(n_-AX,) 



(4.9) 



0<s<t 



- F{Ys-) + F'(n_)AX,]/(AX, < n_) 
+ [F(0) - F{Ys^) + F'(y,_)AX3]/(AX, > Ys-)) 



upon using that AXg < Yg- if and only if Xg < Sg- so that AS'g = 0, and 
AX^ > Yg- if and only if Xg > Sg- so that 5^ = Xg, that is, Yg = 0. Taking Ey 
on both sides of (4.9), where P^^ denotes a probability measure under which 
Yq = y, and applying the compensation formula (see, e.g., [20], page 475) we 
find that 



(4.10) 



EyF{Yt)-F{y) 



ds 



+ 



[F{Yg - x) - F{Yg) + F'{Yg)x]i^{dx) 
[F{0) - F{Yg) + F'{Yg)x]iy{dx 



Ys 



for all y > 0. The applicability of this formula (see, e.g., [15], page 97) follows 
from the facts that \F'{y) \ < Cy and \F"{y) \ < C for all y > with some C > 
so that the mean value theorem yields the existence of ^g^x S iXs — x, Yg) and 
r]g € {0,Ys) such that 



ds 



<E, 



(4.11) 



\F{Yg - x) - F{Yg) + F'{Yg)x\i^{dx 

POO 

+ / |F(0) - F(Yg) + F'{Yg)x\iy{dx) 
Jys 

+ r {\F'{i^g)\Yg + \F'{Yg)\x)-^dx 

JYs X 
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rt 



<CEy 



2(2 — a) a a — 1 



< -(:^7^ + - + f *^/°E,(5i - < oo 

\2[2 — a) a a—lj 2 

since 2 — a € (0, a) and where we also use the scahng property of X. Dividing 
both sides of (4.10) by t, letting t J, and using the dominated convergence 
theorem, we get 



(4.12) 



LyF(y) = nPiy - x) - F{y) + F'{y)x]v{dx) 



roc roo 

+ [F{0)-F{y)] / u{dx)+F'{y) / xu{dx) 

J 11 J V 



which is easily verified to be equal to the right-hand side of (4.3) for all y > 
upon using (2.4). This completes the proof. □ 

3. It will be shown in Section 5 below that the continuation set C in the 
optimal stopping problem (3.2) always contains the interval [0,e) for some 
e > sufficiently small, so that the optimal stopping point from (3.14) 
is always strictly larger than zero. Moreover, we now show that the value 
function V from (3.2) is smooth from the left at whenever D ^ 0. 

Proposition 5 (Smooth fit). If the optimal stopping point from (3.14) 
is finite, then the value function V from (3.2) is differentiable from the left 
at and we have 

(4.13) VL{z.) = G'{z,). 

Proof. To simplify the notation, let us write b in place of z*. Then 
[0,6) CC and fee D so that V{h) = G{h). Hence, {V {h - e) - V {h))/ {-e) > 
{G{h — e) — G(6))/(— e) for all e > sufficiently small, and letting e J, we 
obtain 

(4.14) ^■^.^^y(6-e)-y(6)^^. 

£4-0 -e 

To derive a reverse inequality, note that the stopping time 

(4.15) cj£ = inf {s > : Z^"^ > b} 

is optimal for V{b — e) under P (recall that Z creeps upward). Hence, by the 
mean value theorem we find that 

V{b-e)-V{b) 

>Eie-^-^G{Z'-n)-He-''''GiZ')) 
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(4.16) =^{e-^"^G'{ie){Zlr-Zlj) 

where ie G (^^jS ^Ij for e G (0, 5). 

We claim that (Tg ^ P-a.s. as e J,0. Indeed, setting 

(4.17) p, = inf{s > : (6 - e) V - X^^,^ > b}, 

(4.18) Te = mi{t>0:{b-e)W St-Xt>b} 

we see that cJe < and = t~^(re) for all e > 0. Since t~^{0+) = it is 
therefore sufficient to show that — )• P-a.s. as e ^0. For this, note that 

(4.19) Ts < mf{t >0:{b-e)-Xt>b} = inf{t >0:Xt<-e} =: 7^ 

and 7£ i P-a.s. as e 4, since the starting point of X is regular for (—00, 0). 
Hence, ^ P-a.s. for e 1 as claimed. 

Dividing both sides of (4.16) by —e, letting e 4, 0, and using the dominated 
convergence theorem [upon noting that < + {Z^ — Z^~'^) <b + ed^^ < 
{b + e)e''' and recalling that G'{z) =pzP-^Fs^{z) <2zP-^'iov ah z > so 
that < e-P''-G'{^e)e'''I{St(^,) <b)< 2e^-P+^>- (b + e)P-^e'-P-^^''- = 2(6 + 
e)P"i <2(6 + l)P-i as e|0], we get 

(4.20) limsup^^^^^£l^^<G'(6). 

£4-0 

Combining (4.14) and (4.20), we see that V is differentiable from the left 
at b and that (4.13) holds as claimed. This completes the proof. □ 

4. Returning to the case when [0,z*) C C with z^, & D, recalling the gen- 
eral fact on the killed Dirichlet problem (which suggests that z 1— >■ V{z) = 
Eze~P'^''* G{Zfj^^ ) should solve L^y = pV in [0, z*) due to the strong Markov 
property of Z; see, e.g., [18], pages 130-132), and making use of the facts 
from Propositions 3-5, we can formulate the following free-boundary prob- 
lem for the value function V defined in (3.2) above: 

(4.21) zV'{z)+aLYV{z)-pV{z) = Q for^G[0,z*), 

(4.22) y(z*)=G(z*) (instantaneous stopping), 

(4.23) V' {z^) = G' {z^) (smooth fit), 

(4.24) y(0)=0 (normal reflection), 

where G (0, 00) is the (unknown) boundary point to be found along with V 
on [0,2;*). Whilst the infinitesimal generator Ly in (4.21) can take any of 
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the three forms (4.3)-(4.5) from Proposition 4, it turns out that the Caputo 
form (4.5) is most convenient for the analysis of the problem (4.21)-(4.24) 
to be performed. 

For this reason, let us rewrite (4.21) in the Caputo form as 

(4.25) .F'(.-) + -^|'^-£|L_,,_pF(,) = „ 

for z G (0, 6] and F : [0, 6] -^M. with h G (0, oo) given and fixed. The proof of 
Proposition 6 below shows that the natural solution space for this equation 
is one-dimensional [once -F'(O) is set to 0]. More precisely, let Sh denote the 
class of functions F : [0, 6] ^ M satisfying the following three conditions: 

(4.26) FGCi[0,6]nC2(0,6], 

(4.27) \F'\z) I = 0(z""=^) as z 4 0, 

(4.28) F'(0) = 0. 

Note that F" is assumed to exist (and be continuous) on (0, 6] but may 
be unbounded (locally at zero). Note also that (4.26)-(4.28) imply that 
|-F'(z)| = 0(z"~^) as z 4- 0- For further reference, let us also recall the fol- 
lowing well-known identity (see, e.g., (3.191) in [10], page 333, and (6.2.2) 
in [1], page 258): 

(4.29) r x^~\z - xr~' dx = 

Jo r(^ + u) 

for ^ > and > 0. 

Proposition 6. The equation (4-^5) has a unique solution F in Sb 
satisfying 

(4.30) F(0) = ao, 

whenever cq € M is given and fixed. Moreover, the following explicit repre- 
sentation is valid: 

°^ 1 f r)\ 7°-'"' 

(4.31) F{z) = aQy^ -- r 

^ ' ^' °^^(-cr(-a))"V «y„r(an + l) 

for z G [0, h] where {q)n = q{q + 1) •••((? + n — 1) for n > 1 and {q)o = 1 with 
q = -p/a. 

Proof. 1. Uniqueness. We will establish the uniqueness of solution by 
reducing the integro-differential equation (4.25) to a Volterra integral equa- 
tion of the second kind. For this, let us introduce the following substitution 
in (4.25): 

F"(x) 

(4.32) .(.) = j; (j4i^<ix 
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for z > upon extending F from [0, h] to a bounded function on (0, oo) 
with bounded support in Let L[/](A) = e~^^ f{x) dx denote the 
Laplace transform of a function / : M_|_ — M with A > 0, and let denote 
the inverse Laplace transform. By (4.27) and (4.29) we see that L[99](A) 
is well defined and finite for all A > 0. Applying first L and then on 
both sides of (4.32) using the well-known properties (i) h[f^ fi{y)f2{x — 
y)dy]{X)=Uh]iX)Uh]{X), (ii) L[r](A) = A2L[/](A)-A/(0)-/'(0) and 
(iii) h[xP]{\)=T{p + l)/XP+^ for p> -1, one finds using (4.28) that 

(4.33) F{z) = — -i- Hz - x)"- V(x) dx + F(0) 

r(a)r(2 - a) Jo 

for z £ {0,b]. Inserting this expression back into (4.25), we obtain 

(4.34) [ K{z,x)ip{x)dx + ip{z) ='>p, 
Jo 

where K and ^ are given by 

(4.35) ^(-■^)- wtrTp' /"/"t-."' - 

cr(Q)r(2 — a) [z — xy ° 

(4.36) V = ^^^^F(0) 

c 

for z € (0,6] and x S (0,z). We may now recognise (4.34) as a Volterra in- 
tegral equation of the second kind with a weakly singular kernel K [the 
kernel is said to be weakly singular since the exponent 2 — a in the singular 
term (z — x)^~" belongs to the interval (0,1)]. Moreover, since ^ defines 
a bounded function on [0,6], it is well known (see, e.g., [13], Theorem 7, 
page 35) that the equation (4.34) can have at most one solution 99 (in the 
class of locally integrable functions), and by means of the identity (4.33) this 
fact translates directly into the uniqueness of solution for (4.25) as claimed. 
This completes the first part of the proof. 

2. Existence. Seeking a solution to (4.25) of the form 



(4.37) F{z) = Y,anZ-^ 



n=0 

and inserting it into (4.25) upon differentiating and integrating formally 
term by term and making use of the identity (4.29), a lengthy but straight- 
forward calculation shows that /? = a, 7 = and the series coefficients satisfy 

I fp \ Tian + l) 

(4.38) On+l = —, r n\ —— . . fln 

cT[—a)\a / r(Q(n -I- 1) -I- 1) 

for n = 0, 1, . . . . This yields the candidate series representation (4.31). More- 
over, setting bn = (l/(—cr(— «))"■)(— p/a)n(^""/r(an-|- 1)) for n > 1 and us- 
ing the well-known fact that T{an + 1) /T{a{n + 1) + 1) ~ (an)~" as ?i — )■ 00 
(see (6.1.47) in [1], page 257), it is easily verified that 6„+i/6„ — > as n — )• 00. 
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Hence, by the ratio test we can conclude that the series in (4.31) converges 
absolutely for every z € [0,b]. A direct verification also shows that the func- 
tion F defined by the series in (4.31) belongs to Sb- These facts justify the 
formal steps leading to (4.38) above, and the proof is complete. □ 

5. Before we continue our analysis of the free-boundary problem (4.21)- 
(4.24), let us make precise the following consequence of Ito's formula and 
the optional sampling theorem. Note that G satisfies both (4.39) and (4.40) 
below since |G"(z)| = 0(zP+"-3) as z i and \G" {z)\ = 0{zP-^) as z too. 
This is easily seen upon recalling the expression for G" from the proof of (3.3) 
above and using the asymptotic relations (5.14), (5.15) and (5.17) below. 
Recall also that F from Proposition 6 satisfies (4.39) below. 

Proposition 7. Let F : M+ — M 6e a function from [0, oo) n C^(0, oo) 
satisfying 

(4.39) |F"(z)| =0(z"-2) aszlO and F'(0) = 0. 
Let a be a stopping time of Z such that either 

(4.40) \F"{z)\ = 0{z^) as z-\ oo for some 15 < a -2 

and a <k for some k>l, or a < am for some m > 1 where am = inf{s > : 
Zs = m}. Then the following identity holds: 

(4.41) E,e-P-F(Z,) = F(z) + E, f e-P%l.zF - pF){Zs) ds 

Jo 

for all z G M+ . 

Proof. Under with z € IR+ by Ito's formula, we get 

e-P'F{Z,) = F{z) -p r e-P'FiZr) dr 

(4.42) ^ 

+ / e-P'^F'{Zr)Zrdr + Ms + Js, 
Jo 

where Mg = — e^~P~^^^^'F'{Zr-)dX^(^^.^ is a local martingale and Jg = 
Eo<r<.e-^^-(^(^r) - F{Zr-) + e'^F'(Z,_ ) AXj(,,) ) for s > [upon noting 
that dZr = Zr dr-\-e'' dY^j-^,) from (2.15) above]. Note also that Z is a quadratic 
pure jump semimartingale (i.e., [Z, Z^ = 0) for the reasons outlined follow- 
ing (4.7) above. Note further that similarly to (4.9) we find that 

Js=Y. e-^'([i^(e^>^tM--e'A^*M) 

0<r<s 

- F(e^yi(,)_) + e'^F'(e'-yi(,)_)AXi(,)] 



PREDICTING THE ULTIMATE SUPREMUM 19 

(4.43) X I{AX,^r) <Ytiry) 

+ [F(0) - F(e^y,(,)_) + e^F'(e^yi(,)„)AX,(,)] 

Xl{AXt^,.)>Yt^r)-)) 

upon using that AXi(^r) ^ ^t(r)- if o'^ly ^t(r) ^ Sf(^r)- so that A5f = 0, 
and AX((^) > ^((r)- if and only if X^-^^j.-^ > ^^(r)- so that 5((r) = ^4(^)1 that 
is, yt{r) = 0. Setting u = t{r) this further reads 

0<i'<i(s) 

- F(e*"'(")y„_) + e*"'(^)F'(e*"'(^)y„_)AX„] 

(4.44) 

x/(AX„<y„_) 

+ [F(0) - F(e*~'(^')y„_) + e*''(^')F'(e*~'(^')y„_)AX^] 

x/(Ax,>y,_)). 

The compensator i^' of J is given by 







(4.45) 



- F(e*"'(^')n) + e*"'(^)F'(e*^'(")y„)x]z/(dx) 

/■oo 

+ / [F(0)-F(e*"'(^)y^,) + e*"'(^')F'(e*"'(^)y„)x]i/((ix) . 



Setting r = t~^(v) and y = e'x we see that dv = ae""*" and dx = dy so 
that v{dx) = cdx/x^+°' = {e'^'^+°'> cdx) /y^+°' = {e°"' cdy) /y^+°' = e°V((iy). 
This shows that 

Ks = aj^ e~P' dr l^j^^ [F{Zr - y) - F{Zr) + F' {Zr)y]u{dy) 

(4.46) + ^"[^(0) - F{Zr) + F'{Zr)y]v{dy)^ 

= a r e-P'-LYF{Zr )dr 
Jo 

upon recalling the argument following (4.12) above to obtain the final equal- 
ity [where Lyi^ denotes the action of Ly on F given by the right-hand side 
of (4.3)-(4.5)]. 
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If m > 1 is given and fixed then (4.39) implies tlie existence of C > such 
that \F'{z)\ < Cz""-^ and \F"{z)\ < Cz"^-^ for all z G (0,?n]. This combined 
with the mean value theorem yields the existence of S,r,y £ (Zr — y,Zr) and 
r]r € {0, Zr) such that 

dr\ 1^^ \F{Zr -y)- F{Zr) + F'{Zr)yMdy) 



sAa-n 



-pr 



+ 



POO 

\F{0)-F{Zr) + F'{Zr)y\Hdy) 



sACn 



-pr 



dr 



(4.47) 



+ 



\F'{r]r)\Zr + \F'{Zr)\y)- 



,1+a 



dy 



<cE, 



sAcr„ 



-pr , 



c 



dr{- (Zr-y) 



a-2„.l-a 







2 

+ cz:, 



dy 



-l-a 



dy + CZ^ 



a-l 



Zr 



y-'^dy 



:c(^r(2-a)r(a-l) + ^ + ^ 
2 a a — I 



SACTry 



e-P'' dr 



< oo 



upon using (4.29) in the final equality. It follows that A^sact^ •= JsAa^n ~ 
KsA<Tm, is a martingale under for s > (see, e.g., [15], page 97). This 
shows that := J — X is a local martingale [with {<Jm)m>i as a localization 
sequence of stopping times] . 

Let o" be a stopping time of Z such that a < am for some m > 1. Choose 
a localization sequence of stopping times {pn)n>i for the local martingale M. 
Subtracting and adding Kg on the right-hand side of (4.42), replacing s by 
a Apn, taking on both sides and applying the optional sampling theorem, 
we obtain 

(4.48) E,e-P('^^^")F(Z,ApJ = F{z) + / e-P'' {hzF - pF){Zr) dr 

Jo 

for all z € M+ and all n > 1 [upon recalling (4.46) and the action of 
in (4.1) above]. Moreover, it is easily seen from (4.5) using (4.39) and (4.29) 
that z i-7> LyF(z) is bounded on [0,m] (and so are F and F' by continuity). 
Letting n ^ oo in (4.48) and using the dominated convergence theorem we 
see that (4.41) holds as claimed in this case. 

Let us now assume that (4.40) holds with a <k for some k>l. Choose 
again a localization sequence of stopping times {pn)n>i, however, this time 
for both the local martingale M and and the local martingale N. Subtracting 
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and adding Kg on the right-hand side of (4.42), replacing s by a A pn, ta- 
king Ez on both sides and applying the optional sampling theorem, we again 
obtain (4.48) for all z € M+ and all n > 1. Moreover, it is easily seen from 
(4.5) using (4.39)+(4.40) and (4.29) that |LyF(z)| < C3(l + 0^+^"") for all 
z € M+ with some C3 > 0. Likewise, it is easily verified that (4.39) and (4.40) 
imply that \F{z)\ < C^il + 2^+^) and \F'{z)\ < C^il + z^^^) for all zeR+ 
with some C4 > and C5 > 0. Hence, we see that there exists Cq > such 
that 

\F{ZI)\ + \{hzF-pF){ZI)\ 

(4.49) <Ce{l + {Z'sf^^) 

<C6(l + e'=(^+2)(z + 5i-/i)'^+2) 

for all s € [0,k] where the right-hand side defines an integrable random vari- 
able since /3 + 2 G (0, a) . (Note that without loss of generality, we can assume 
that 13 is close enough to a — 2 so that /3 + 2 > 0.) Letting n — >■ 00 in (4.48) 
and using the dominated convergence theorem (twice) we see that (4.41) 
holds as claimed. This completes the proof. □ 

6. We now establish a remarkable probabilistic representation of the global 
solution (4.31) to the equation (4.25). For this, let us set 

(4.50) Vi{z) = E{zV Si-Xi)P 

for all zeR+. From (3.13), we see formally that Vi{z) = E^e-P'^°°G(Z^^) for 
all z G M+ where Uoo = inf{s >0:Zs = 00}, and this suggests that z 1-^ Viiz) 
should solve the equation (4.25) on R_|_. This can be derived rigourously as 
follows. 



Proposition 8. Let Fi denote the global solution (4-31 ) to (4-25) on \ 
with -Fi(O) = 1. Then the following identity holds: 

(4.51) Vi{z) = aiFi{z) 

for all z G where the constant oi is given explicitly by 

(4.52) a, = a{cr{-a)r/^^^ 

T{p/a) 



Proof. 1. We first show that the identity (4.51) holds with some cons- 
tant ai>0. For this, fix an arbitrary zi>,0, set F[z) = aFi{z) for 2;GM+ 
where a = Vi{zi) / Fi{zi) , and consider = inf{s >Q:Zs = zi}. Then by 
(4.41) and (4.25), we find that 



(4.53) F{z) = E,e-f^^iF(Z,,^ ) = F(zi)E,e-?"^^i = yi(zi)E,e-P'^^i 
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for all z G [0, z-\\. In addition, consider On = inf{s > : = n} and set 

(4.54) y"(z) = E,e"P'^"G(Z,J 

for n> zi and z G [0, zi]. Note that (3.13) implies that V"{z) — > Vi{z) as 
n — >■ (X) for all z G [0, zi]. Fixing n > zi and applying the strong Markov 
property of Z at we find that 

(4.55) = E,(e-P'^"°^-i-^"^^i+f^^iG(Z.J o 0,^.^ |7-,(,^^)) 

= e^'-iE,(e-*'-G(Z.J|J-,(,^^)) 

for all z G [0,2;i]. Multiplying both sides by e"^'^''! and then taking E^, we 
get 

(4.56) F"(zi)E,e-f"^i =y"(z) 

for all z G [0, zi] and n> zi. Letting n — ?• oo we obtain 

(4.57) yi(zi)E,e-P'^^-i =yi(z) 

for all z G [0,zi]. Comparing (4.57) with (4.53), we see that Vi{z) = F{z) for 
all z G [0,2;i]. Since zi > was arbitrary this establishes (4.51) with some 
constant oi > 0. 

2. To derive (4.52), we may apply the Laplace transform L on both sides 
of (4.25) where F{z) = Vi{z) = aiFi{z) for z G M+ so that ai = 1^(0). Us- 
ing the well-known properties (i)-(iii) recalled following (4.32) above and 
(iv) L[zF'(z)](A) = -AL[F]'(A) - L[F](A) for A > 0, it can be verified using 
(4.28) that this leads to 

Mi^r(A) + (^-^^^A"-)L[F](A) 

(4.58) 

= -F(0)^^^^^A""2 
a — 1 

for A > 0. Solving this equation under L[F](A) — )• as A ^ oo [this condition 
is satisfied since F{z) = Vi{z) ~ z^ as z — >• oo by (4.50) above] we find that 

(4-59) L[F](A) = r(l+p/a,cr(-a)A-) 

for A > 0, where V{a,x) = y'^~^e~y dy denotes the incomplete gamma 
function for a > and x > 0. Since z i— ?> F{z) is increasing [by (4.50) above], 
we can use the Tauberian monotone density theorem (see, e.g., [15], Theo- 
rem 5.14, page 127) which states that (i) L[F](A) ~ 1\~p as A J, if and only 
if (n) F{z) ~ {i/V{p))zP~^ as z t oo where p > and £ > 0. From (4.59), we 
see that (i) is satisfied with p = 1 +p and £ = (F(0)/(cr(-a))P/°)r(l +p/a) 
so that (ii) yields (4.52) since F[z) = Vi{z) ~ z^ as z — >• oo. This completes 
the proof. □ 
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5. Predicting the ultimate supremum. 1. We will begin by connecting 
our findings on the free-boundary problem from the previous section to the 
value function from (3.2). 

Proposition 9. // the optimal stopping point from (3.14) is finite, 
then the value function V from (3.2) coincides on [0,z*] with F from (4-31) 
where oq is set to ^(0). In terms of the function V\ from (4-50), this reads 
as follows: 

(5.1) V{z)=aVi{z) 

for all z € [0,z*] where a = V{0)/ai G (0, 1) and ai is given by (4-52) above. 
If the optimal stopping point is not finite (i.e., the optimal stopping set D 
is empty), then 

(5.2) V{z) = Vi{z) 
for all z G M+ . 

Proof. If < oo then 

(5.3) V{z) = E,e-P''^*G{Z^J = y(z,)E,e-f'^- 

for all z € [0, z*]. Moreover, if we set F{z) = aQFi{z) for all z G M_|. with 
ao = V{z^:)/Fi{z^) then by (4.41) and (4.25), we have 

F{z) = E,e-P'^-F(Z,_J = F(z,)E,e-?''^- 

(5.4) 

= y(z.)E,e-f'^- 

for all z G [0,z*]. Comparing (5.3) and (5.4), we see that V{z) = F[z) for all 
z G [0,z*]. Hence oq = ^(0) and this establishes (5.1) upon recalling (4.51). 
If z* = oo then (5.2) follows from (3.13) above. This completes the proof. 
□ 

From (5.1) and (5.2), we see that the value function V \s a constant multi- 
ple of the function Vi from (4.50) up to the first contact point with G (when 
starting from and moving toward oo in the state space). The unknown 
constant needs to be chosen so that the contact with G occurs smoothly. 
Since V <Vi this leads the following criterion for D to be nonempty: 

(5.5) < oo if and only if 3zi G M+ such that ^1(2:1) > G{zi) 
or equivalently, the following criterion for D to be empty: 

(5.6) 2* = 00 if and only if Vi{z) < G{z) for all z G M+. 
We will continue our analysis by examining when (5.5) holds. 
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2. Consider the function H : [0, oo) M defined by 

(5.7) H{z) = {hzG-pG){z) 

for z > where H{0) := H{0+) exists by (5.11) below. Recall that (4.41) 
reads 

(5.8) E^e-f"G(Z<,) = G{z) + E, T e~P'H{Zs) ds 

Jo 

for z € where a is any stopping time of Z like in Proposition 7. Set 

(5.9) N = {ze[0,oo):H{z)<0} and P = {z e [0,oo) : H{z) >0}. 
Then the fohowing two inclusions are vahd: 

(5.10) NQG and D <Z P. 

Indeed, to show the first inckision (the second one then being obvious) 
take any z & N and choose e > smah enough such that {z — e,z + e)ri M+ C 
A'^ (note that N is open in M+). Inserting the stopping time a,. = inf {s > 
0:Zs ^ {z - e,z + e)} into (5.8), we see that E^e~P'^=G(Z^J < G{z) since 
H{Zs) < for s S [0,(Te). Hence, z belongs to G as claimed. 

3. Motivated by the important role that the function H plays in the 
optimal stopping problem (3.2), we now determine its asymptotic behav- 
ior at zero and infinity. Note that (5.11) below and (5.10) above imply 
(since H is continuous) that the continuation set G always contains the 
interval [0, e) for some e > sufficiently small so that the optimal stopping 
point from (3.14) is always strictly larger than zero. 

Proposition 10. The following relations are valid: 

(5.11) \imH{z) = -pG{0) = -pESf<0, 

(5.12) hm z--PH{z) = —^^—-(rip -a)- r(p)r(l - a)), 
ztoo L [p — a + 1) 

Proof. Since G'{z) = pz^'^Fs^iz) and G"{z) = p{p - 1)zP-^Fs^{z) + 
pzP~^ fsi{z) we see by (4.1) and (4.5) that 

Hiz) = zG'iz) + ^J^'^^^d.-pGiz) 
= pzPFs,{z) 

(5.13) 

c r p{p - l)xP-^Fs, (x) + pxP-^fs, jx) 
a- I Jo (z- x)""-^ 

-pG{z) 



PREDICTING THE ULTIMATE SUPREMUM 25 

for z > 0. Recall that the following asymptotic relations are valid (see [2], 
Corollary 3): 

^'■''^ ^"-^^^^ (cr(-a))i-v.r(a-i)r(i/a) 
^"-^^^^ (cr(-a))i-v»r(a)r(i/a) 

Using (5.14) and (5.15) together with (4.29) it is readily verified that the 
integral in (5.13) tends to as z],0. This easily yields the first equality 
in (5.11) and the second equality follows from (3.3). 

Moreover, using (2.12) above we can further rewrite (5.13) as follows: 

H{z) = -pzP{l-FsAz)) 

cp(p - 1) xP""^ 

+ ^ ' / 7 ^ — jdx 

a -I Jo (z-x)°-i 

POO 

(5.16) -p il-Fs,{x'^n)dx 

JzP 

a-l Jo (z-x)"-^ 

a-l Jo (z-x)"-i 

for z > 0. Recall that the following asymptotic relations are valid (cf. [2,5, 
16]): 

(5.17) fsM 

(5.18) i-F5,(z)~^ as z too. 

az" 

Using (5.17) and (5.18) together with (4.29) it is somewhat lengthy but still 
straightforward to verify that the final two integrals in (5.16) are o(zP~°') as 
z oo, whilst the first three terms in (5.16) multiplied by z"~p converge to 
the constant on the right-hand side of (5.12) as z — > oo. This completes the 
proof. □ 

4. Motivated by the identity (5.12) let us consider the function i defined 

by 

(5.19) iia,p) = ^,^-^—^{r{p -a)- r(p)r(l - a)) 

for a € (1,2) and p G (l,a). A direct examination of the right-hand side 
in (5.19) shows that there exist a* € (1,2) (equal to 1.57 approximately) 
and a strictly increasing function : (a*, 2) — )• (1, 2) satisfying p^{a*+) = 1, 
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p*(2— ) = 2 andp*(a) < a for a G (a*, 2) such that (i) i{a,p) > if a S (a*, 2) 
and p G (l,p*(a)); (ii) £{a,p) < if either a € (1, a*) and p G (1, a) or q G 
[a^,,2) and p G (p*(a),Q); and (iii) £(a,p*(a)) = for a G (a^=,2). Note that 
the properties (i)-(iii) do not depend on the value of the constant c in (2.4). 
Recall also from (5.12) above that 

(5.20) i{a,p) = lmiz''-PH{z) 

ztoo 

for all a G (1,2) and p G (l,a). In view of (5.10) this suggests that the sign 
of £ plays an important role in the problem (3.2). 

Building on the facts presented in the previous sections, and extending 
these arguments further in the proof below, we can now present the main 
result of the paper. It should be recalled in the statement below that the 
function Vi can be expressed probabilistically by (4.50) and analytically by 
(4.51)+(4.52) [where Fi is given by (4.31) with oq = 1], and the probabilistic 
and analytic representations of the function G are given in (2.12) above 
(upon recalling that Fs-^ admits an explicit series representation as shown 
in [2], Theorem 1). 

Theorem 11. I. If a € (a*, 2) and p G (l,p*(a)) then there exists G 
(0,oo) such that the stopping time (3.14) is optimal in the problem (3.2) 
under for z G [0,z*]. The optimal stopping point can he characterized 
as the minimal z G (0, cxd) for which 

(5.21) P*Vi{z)U=,^=G{z)U=,„ 

where /3* G (0,1) is the minimal j3 G (0,1) for which (5.21) has at least one 
root z G (0, oo) . The optimal z^ and satisfy the smooth fit condition 

(5.22) l3,Vl{z)U=,,=G'{z)U=,^. 

The value function from (3.2) is given by V{z) = /3*Vi(z) = l3^E{z\/ Si — Xi^ 
for z G [0, z*] . 

II. The stopping time (1.3) is optimal in the problem (2.3) and the value 
from (2.3) is given byV = TP/''l3^Vi{0)=TP/''(3^E{Si-Xi)P = TP/''f3^a{cx 
T{-a))P/T{p)/r{p/a). 

Proof. Since part II follows from part I as discussed in Sections 2 
and 3 above, it is enough to prove part I. For this, we will first show that 
the assumptions a G (a*, 2) and p G (l,p*(a)) imply the existence of zi > 
(large enough) such that 

(5.23) Vi{z)>G{z) 

for all z > zi. We will then show how the knowledge of (5.23) combined with 
the properties and facts about Vi and G derived in the previous sections yield 
the existence of /3* and z* satisfying the remaining statements of part I. 
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1. To prove (5.23), recall that the identity (4.41) is applicable to G in 
place of F with a = n for n > 1. Letting n ^ oo in this identity, using (3.13) 
combined with the fact that each e~P"'G{Z^) is dominated by {z V Si — 
IiY + ESf which clearly has finite expectation, as well as the fact that the 
function H is bounded (by the result of Proposition 10), it follows by the 
dominated convergence theorem that 

/>oo 

(5.24) E(zV5i-Xi)P = G(2) + E / e~P'H{Zl)ds 



for all z > 0. Recognizing the left-hand side of (5.24) as Vi(z), we see 
that (5.23) will be established if we show the existence of zi > (large 
enough) such that 

poo 

(5.25) I{z) := E / e-P'H{Zl) ds > 

Jo 

for all z> zi. 

To show (5.25) recall from (i) following (5.19) above that i := i{a,p) 
in (5.20) is strictly positive when a € (a*, 2) and p G (l,p*(a)) are given and 
fixed. Hence for any given and fixed e > (small) there exists z^ > (large) 
such that 

(5.26) z''-PH{z)>i-e 
for all z> Ze- Consider 

/•oo 

(5.27) J{z) := E / e-P'H{Zl)I{Zl < z^) ds, 

Jo 

/•oo 

(5.28) K{z) := E / e-P'H{Zl)I{Zl > Ze) ds 

Jo 

and note that I{z) = J{z) + K{z) for all z > 0. 

Let M > be large enough so that |-ff(z)| < M for all z > 0. Then we 
have 

/•oo 

(5.29) \J{z)\<M e-P'P{Zl <ze)ds 

Jo 

for all z > 0. Moreover, by (5.18) we see that 

P{Zl < Ze) = P{e%Z V %) - X,(,)) < Ze) 

< P(Z V 5t(s) - S't(^) < Ze) 

(5.30) < P(z - 5i(,) < Ze) 

< P{Si >Z-Ze) 
<N-{Z-Ze)-'' 

a 
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for all z > Ze with some > large enough. Combining (5.29) and (5.30) 
we find that 

(5.31) |j(^)|<^(^_^^)-" 

pa 

for all z> z^. 

On the other hand, by (5.26) we see that 

K{z) = E / e-P'H{ZI)I{ZI > Ze) ds 
Jo 

j-OO 

>{i-e) e-P^E[{Z!r~'^I{Z^,>z,)]ds 
Jo 

POO 

= {e-e) e-°^E[(zV%)-Xi(,))P-"/(Z|>z,)](is 

(5.32) ^° 



>U- e)zP-'' 



ds 



poo 

/ e-'^'EKlV Si-Ii)P^'^IiZI>Ze)]ds 
Jo 



> ^— ^(E(l V 5i - hf-" - S)zP-" 
a 

for all z>lVzs, where in the second last inequality we use that 

(5.33) lv:^-^<lV^-:^<lV5i-/i 

z z z z 

for all s > and z>l, and in the last inequality we use that 

/•oo 

lim / e-"^E[(lV5i-/i)P-"/(Z|>z,)]ds 

(5.34) 

= -E(lV5i-/i)P-"<oo 
a 

by the dominated convergence theorem since Z| — ?> oo as z ^ oo [from (5.34) 
we see that for given 6 € (0, E(l V 5i — /i)^"") there exists zs > such that 
the final inequality in (5.32) holds for all z > zs]. Since the right-hand side 
in (5.31) tends faster to zero than the right-hand side in (5.32) as z oo, we 
see that (5.23) holds with some zi > large enough as claimed. 

2. We now establish the existence of /3* and satisfying the remaining 
statements of part I. For this, recall that (5.23) holds for z = zi so that 
for some /3i G (0,1) sufficiently close to 1 we have f3iVi{zi) > G{zi). Since 
/3iVi{z) ~ /3izP < z^' ~ G{z) as z ^ oo we also see that there exists Z2 > z\ 
such that /3iVi(z) < G(^z') for all z > 22- This shows that for some /3o G 
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(0, 1) sufficiently close to we have /3oVi(z) < G{z) for all z > [recall that 
Vi{0) = E(5i — Xi)P > and that Vi is increasing]. It follows therefore by 
continuity that there exists the smallest /3* € (/3o,/3i) C (0,1) such that the 
set A = {z e M+|/3*Fl(z) = G{z)} is nonempty so that l3Vi{z) < G{z) for 
all z G if /3 G (0,/3*). Setting w^, = miA we see that w^: belongs to A by 
continuity so that (5.21) holds for z = w*. Moreover, since Vi{zi) > G{zi) we 
know by (5.5) that z^ = inf D < oo so that by (5.1) we have V{z) = a^Vi{z) 
for all z G [0,z*] with some a* G (0, 1). By the construction of /3* and it 
follows therefore that /3* < a* and w^> z^. If either < a* or equivalently 
> z*, then since P*Vi{z) = Ee"?"^™* (/3*Vi)(Z|^ J for all z G [0,w*] by the 
result of Proposition 7, and this further equals Ee~^'^™* G{Z^^ ) for all 2; G 
[0,tt;*] by definition of aw,, we see that /3*Vi(0) > V{0) while at the same 
time /3^Vi{0) < a*Vi(0) = V^(0) which is a contradiction. Thus /3^, = a* and 
If* = 2:* so that V{z) = f3^Vi{z) for all z G [0,z*] as claimed. The smooth fit 
condition (5.22) then follows by the result of Proposition 5. This completes 
the proof of part I whence part II follows as discussed above. □ 

5. In the final part of this section, we briefly consider the case when the 
hypotheses of Theorem 11 are not satisfied. 

Proposition 12. If either a G (l,a*) orpG (p*(a),a), then there exists 
zi> large enough such that Vi{z) < G{z) for all z>zi. 

Proof. This can be proved in exactly the same way as (5.23) above 
upon noting that I :=l{a,p) in (5.20) is strictly negative when either a G 
(1, a*) or p G {p^{a),a) and replacing (5.26) with z'^~^H{z) < £+e for all z > 
Ze- This leads to (5.31) without changes and (5.32) holds with the inequalities 
reversed since £ + e < in this case. Different rates of convergence in the 
resulting inequalities then complete the proof just as above. □ 

It follows from the result of Proposition 12 that the continuation set C 
contains the interval [2:1,00) for some 21 > large enough when either a G 
(l,a,,) or p G (p*(a),Q!). It shows that the stopping time (1.3) can no longer 
be optimal in this case (in the sense that it is not optimal to stop at t G [0, T) 
when St — Xt is sufficiently large). This stands in sharp contrast with the 
Brownian motion case (formally corresponding to a = 2) where it is optimal 
to stop in such a case. Recall also that the continuation set C always contains 
the interval [0, e) for some e > sufficiently small so that the stopping set D 
must be contained in [s,zi — 6] for some 5 > 0. We do not know whether 
Vi (2) < G{z) holds for all 2 G M+ in this case, or equivalently, whether the 
stopping set D is empty [recall (5.6) above]. This is an interesting open 
question. We refer to [6], Figure 1, for a related phenomenon in the presence 
of strictly positive drifts and the absence of jumps. 
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APPENDIX 

In this section, we determine the action of the infinitesimal generator of 
the reflected process Y = S — X when X is a general (strictly) stable Levy 
process (see [25]). Set 

(A.l) Uo.{dx) = >Q)dx+ ^J~y^J {x < 0) dx, 

where c+ and c_ are nonnegative constants (not both zero) and a € (0,2). 
For a = 1, the two constants need to be identical (see, e.g., [21], pages 86 
and 87), so that 

(A.2) ui (dx) = -^/(x / 0) dx 

x^ 

with c > 0. Recall that C^(M_|_) denotes the class of twice continuously dif- 
ferentiable functions F : R_|_ — >■ M such that F' and F" are bounded on R_|_ . 

Proposition 13. Let X = {Xt)t>o be a stable Levy process of index 
a € (1,2) whose characteristic function is given by 

Ee'^^' = exp(^t j (e^^^ - 1 - iXx)u^{dx)^ 

(A.3) 

^ ^(c+{-i\r+c_{ixr)r{-a)t 

for A € M and t >0. Then the infinitesimal generator Ly of the reflected 
process Y = S — X takes any of the following three forms for y > given 
and fixed: 

ltd 's form 

ry 

LYF{y)= {F{y-x)-F{y)+F'{y)x)-^dx 
Jo X 

(^^4) ^ c+{F{0)-F{y)) ^ c+F'{y) 



ay"' (a - l)y°~^ 

+ / (F(y + x) - Fiy) - F'{y)x)^dx, 







x^ 



Riemann-Liouville 's form 



. ^ C+ d^ fy F(x) , C+F(0) 

^YF{y) = —^—. ^->dx+ + 



aia — 1) dy"^ Jn (y — x)" ^ ay" 
(A.5) 

c- /-°° F(x) 

^a{a-l)dy^Jy {x - y)-^ 
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Caputo 's form 
A.6 

+ —, r / V^- dx, 

a{a-l)Jy 

whenever F € C^(M+) satisfies 

(A. 7) -F'(0+) = (normal reflection) 

with \F"{y) I = 0{y"') as y — )• oo /or some 7 < a — 2 /as weZZ as \F{y) \ = 0{y^) 
as y ^ 00 for some 5 < a — 2 in (A. 5) above]. 

Proof. As in the proof of Proposition 4, it is enough to derive (A. 4). 
For this, fix t > and note that by Ito's formula we have 

F{Yt) = F{Yo) + [' F'{Ys-){dSs-dXs) 

(A.8) 

+ ^ {F{Ys) - FiY,_) - F'iY,_)iASs - AX,)) 

0<s<t 

since [1^,5^]'^ = for the same reasons as in (4.7). Letting Ss = Sg + Sg be 
the decomposition of s Sg into continuous and discontinuous parts, and 
noting that dSg = ASg, we see that (A.8) simphfies to 

F{Yt) = F{Yo) + Mt+ f F'{Yg^)dS'g 

(A.9) 

+ {F{Yg^+AYg)-F{Yg^) + F'{Yg^)AXg), 

0<s<t 

where Mt = - £ F'{Yg_) dXg. Since F' is bounded the same argument as 
following (4.8) above shows that M is a martingale. If s belongs to the 
support of dSg in [0,t], then either Sg_^ < Sg and therefore < 5^ for 
e > implying Yg^ = 0, or Sg < Sg^^ and therefore Sg < Sg+e for e > 
implying Yg = 0. Since there could be at most countably many s in [0,t] 
for which Yg 7^ 1^-, it follows using (A. 7) that the integral with respect 
to dSg in (A.9) is zero. Moreover, the right-hand side of (A.9) can further 
be rewritten as follows: 

F{Yt) = F{Yo) + Mt 

+ J] mYg.-AXg)-F{Yg_) 

0<s<t 

(A.IO) 

+ F'{Yg^)AXg]I{AXg<Yg^) 
+ [F(0) - F{Yg_) + F'{Yg_)AXg]I{AXg > Yg_)) 
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using the same arguments as in (4.9) above. Taking Ey on both sides of (A. 10), 
where Py denotes a probabiUty measure under which Yq = y, and applying 
the compensation formula (see, e.g., [20], page 475) we find that 
EyF{Yt)-F{y) 



(A.ll) 



ds 



[F{Ys -x)- F{Ys) + F'{Ys)x]i^o.idx) 

) 

POO 

+ / [F{0) - F{Ys) + F'{Y,)xMdx: 
JYs 



for all y > 0. The applicability of this formula (see, e.g., [15], page 97) follows 
from the facts that \F'{y)\ < C and \F"{y)\ < C for all y > with some C > 
so that the mean value theorem yields the existence of ^s,x ^ {Yg, Yg + x) and 
ils,x £ (ys,Ys + x) such that 



[ ds( [ \F{Ys-x)-F{Ys)+F'{Ys)x\ua{dx 

Jo \J~OD 

j^dsl^j^ \F{Yg + x)-F{Ys)-F'{YM-l 



(A.12) 



= E 
< E 



■ dx 



ds 



dx 



\F"{^s,.)\x' 



+ 



\F\r]s^^)\x + \F'{Ys)\x) 



dx 



<cE„ 



ds 



C 



+ 



2C_ 



c(7-3a)C 
2(2-a)(a-l) 



t <oo, 



2(2 -a) 

where the remaining two integrals (from to 1^ and from 1^ to cxo) can be 
controlled (bound from above) in exactly the same way as in (4.11) above. 
Dividing both sides of (A.ll) by t, letting t J, and using the dominated 
convergence theorem, we get 

LyF(y)= f [F{y-x)-F{y)+F'{y)x\vo.{dx) 

(A.13) 



+ [F(0) - Fiy)] / u^dx) + F'{y) / xz/„(dx), 

J V J V 



which is easily verified to be equal to the right-hand side of (A. 4) for all 
y > upon using (A.l). This completes the proof. □ 

Proposition 14. Let X = {Xt)t>Q he a stable Levy process of index 
a € (0, 1) whose characteristic function is given by 



Ee 



iXXt 



(A.14) 



exp( t / (e*^^ - l)iya{dx) 



^{c+(-jA)"+c_{jA)")r(-a)t 
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for A S M and t >0. Then the infinitesimal generator Ly of the reflected 
process Y = S — X takes any of the following three forms for y > given 
and fixed: 

ltd 's form 

c+ , c+(F(0)-F(y)) 



LyF(y)= [\F{y-x)-F{y))^dx + 
In X 

(A.15) '° 

^ ' poo 

+ / {F{y + x)-F 
Jo 

Riemann-Liouville 's form 



ay"" 



^ , c+ d fy Fix) , c+F(0) 

(A.16) 



ay" 



H ^ / 7 — —^dx, 

a dy J (x-y)" 



(A.17) LyF(y) = dx+^-f 



■ dx. 



Caputo 's form 

c+ fy F'{x) , c_ p_F^ 

(x - yY 

whenever F € C^(M+) satisfies \F'{y) \ = 0{y'^) as y ^ oo for some 7 < q — 1 
fas well as \F{y)\ = 0{y^) as y — )> 00 for some S <a — 1 in (A.16) above]. 

Proof. As in the proof of Proposition 4 it is enough to derive (A.15). 
For this, fix t > and note that since X is a pure jump semimartingale with 
bounded variation, we have dXg = AX^ and dSg = AS'g for < s < i, so that 
Ito's formula yields 

(A.18) F{Yt) = F{Y,) + {F{Ys) - F{Y,.)). 

0<s<t 

Proceeding as in (A. 10), taking Ey on both sides of the resulting identity 
and applying the compensation formula (see, e.g., [20], page 475), we find 
that 



EyF{Yt)-F{y) 

f ds( r [F{Ys-x)-F{Ys)]vMx) 

Jo \J-oo 

/•CO 

+ / [F{0) - F{Y,)Mdx 
Jy, 



-y 

rt / rYs 

(A.19) = E, 



for all y > 0. The applicability of this formula (see, e.g., [15], page 97) follows 
from the facts that \F{y) \ < C and |-F'(y)| < C for all y > with some C > 
so that the mean value theorem yields the existence of 2, G (Kj , Kj + x) , 
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^s,x ^ i^s - X, Ys) and r]s G (0, Y^) such that 
rt / rY, 
10 



-y 



(A.20) 



ds 



\F{Ys - x) - F{Ys)\u^idx) 

oo 

+ / \F{0) - F{Y,)\i^Udx) 

f ds ( r \F{Ys + x)- F{Ys)\ dx 
Jo \Jo X ^ 



+ / \FiY,-x)-FiY,)\-j^dx 



< c 



+ /j^(0)-W)l-^o 
j/sl^j\F'{es,.)\x^dx 

+ / \F{Ys + x)-F{Ys)\^dx 
+ £\F'ifs,.)\x^dx 

c 

+ / \F{Y,-x)-F{Y,)\-—^dx 
+ f'\F'{r,s)\Ys-^dxI{Ys<l) 

JYs X 

+ f\m-F{Y,)\^dx 



1 — a a a 



2C 7C\ 

H t < oo 



1 — a a 



upon using that 1 — a € (0, 1) in the final inequahty. Dividing both sides 
of (A. 19) by t, letting t^O and using the dominated convergence theorem, 
we get 



(A.21) 



LyF(y) = r [F{y - x) - F{y)]v^{dx) 

J —oo 



+ [F{0)-F{y)] Uo.{dx 
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which is easily verified to be equal to the right-hand side of (A. 15) for all 
y > upon using (A.l). This completes the proof. □ 

Proposition 15. Let X = {Xt)t>o be a stable Levy process of index 1 
whose characteristic function is given by 

(A.22) Ee^^^* = exp(^t (e^^^ - 1 - iXxL{\x\ < l))ui{dx)^ = e^^l^l^* 

for A € R and t >0. Then the infinitesimal generator Ly of the reflected 
process Y = S — X takes any of the following three forms for y > given 
and fixed: 



y c 



ltd 's form 

hYF{y)= I {F{y-x)-F{y)+F'{y)x)^dx 
Jo ^ 

^ c{F{0)-F{y)) 
V 

(A.23) 

+ / iF{y + x)- F{y) - F'{y)x)-^dx 
Jo 

f°° c 
+ {F{y + x)-F{y))—dx, 

Jy 

Riemann-Liouville 's form 

(i^ p / 1 \ cF(0) 



d^ / 1 \ 

(A.24) Ly = c—^ F{x) log [j^^) dx + 

Caputo 's form 

(A.25) hYF{y) = c[ F" {x)log( -^—-] dx, 

Jo \\y-x\J 

whenever F G C^(M+) satisfies 

(A. 26) F'(0+) = (normal reflection) 

with \F"{y) \ = 0{y"') as y ^ oo for some 7 < — 1 [as well as \F{y) \ = 0{y^) 
as y ^ 00 for some 5 < —1 in (A.24) above]. 

Proof. As in the proof of Proposition 4 it is enough to derive (A. 23). 
Using the same arguments as in (A. 8) and (A. 9), we find that 

F{Yt) = F{Yo)- f F'{Y,.)dXs 

(A.27) 

+ {F{Ys-+^Ys)-F{Ys.) + F'{Ys.)^Xs), 

0<s<t 
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where (Jq F'{Ys^) dXs)t>o is a local martingale. We can no longer claim that 
this process is a martingale, however, we note from (A. 22) that Xt = Mt + At 
with 

(A.28) Ee^^^* = exp {t (e^^^ - 1 - i\x)vi [dx]^ , 

(A.29) Ee*^^' = exp( t / (e*^^' - l)vi{dx)\ 

V J\x\>l J 

from where we see that the (Levy) process M = (Mt)t>o is a martingale 
(whose Levy measure has bounded support) and the bounded variation 
(Levy) process A = {At)t>o is given by 

(A.30) At= AXsI{\AXs\ > 1) 

0<s<t 

for t > 0. From (A.27)-(A.30), we see that 

F{Yt) = F{Yo) - f F'{Ys-) dMs 
Jo 

+ (F(n„ + Ay,)-F(n„) 

0<s<i 

+ F'(y,_)AX,/(|AX,| <1)) 

= F(Yo)- [ F'{Ys^)dMs 
Jo 

+ J2 mYs^-AXs)-FiY,^) 



(A.31) 



0<s<t 

+ F\Y,.)AXsI{\AXs\ < l)]I{AXs < Y,_) 
+ [F(0)-F(n_) 

+ F'{Ys.)AXJ{\AXs\ < l)]I{AXs > Ys^)) 

using the same arguments as in (4.9) above. Since F' is bounded and the 
Levy measure of M has bounded support (implying Esupo<s<t|Afs|'J < oo 
and hence E[M, M]'^/^ < oo for all g > by the BDG inequality) it also fol- 
lows by the BDG inequality (with q = l) that (/q F'{Ys^) dXs)t>o is a mar- 
tingale. Taking Ey on both sides of (A.31), where Py denotes a probability 
measure under which Yo = y, and applying the compensation formula (see, 
e.g., [20], page 475) we find that 

EyF{Yt)-Fiy) 

(A.32) =Ey\ [ ds( [ ^ [F{Ys-x)-F{Ys) + F'{Ys)xI{\x\<l)]ui{dx) 
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/•oo 

+ / [FiO)-F{Y,) + F'{Ys)xI{\x\<l)]u,idx] 
Jys 

for all y > 0. The applicability of this formula (see, e.g., [15], page 97) follows 
from the facts that \F{y)\ <C{l + y^^^), \F' {y)\ <C{y ^y'<+^) and \F"{y)\ < 
C(l A y'*') for all y > with some C > so that the mean value theorem yields 
the existence of ^ G {Ys,Ys + x), G {Ys — x,Ys)) and rjs G {Q^Yg) such 
that 



f ds( f ^ \F{Ys-x)-F[Ys) + F'{Ys)xI{\x\<l)\ui{dx) 

Jo \J-oo 

/•oo 

+ / |F(0) - F{Y,) + F'iYs)xI{\x\ < l)\Mdx) 
Jys 

J^dsi^J^ \FiYs + x)-FiYs)\uiidx) 

+ [ \F{Ys + x)-F{Ys)-F'{Y,)x\iyiidx) 
Jo 

+ / \F{Ys -x)- F{Ys) + F\YM^i{dx)I{Ys < 1) 
Jo 

+ \F{Ys -x)- F{Ys) + F'{Y,)x\u,idx)IiYs > 1) 
Jo 

+ \F{Ys -x)- F{Ys)\vi{dx)I{Ys > 1) 

+ / \F{{))-F{Ys) + F'{Ys)x\yi{dx)I{Ys<l) 

JYs 

\F{{})-F{Y,)\vi{dx) 

ds [2c + {Ys + xy+^)-^dx 

+ ri\F"{esJ\x'^dx + 2 ['hF"{ClJ\x'^dx 
Jo 2 ■ x^ Jq 2 ' x^ 

+ 2C(1 + y7+2) / ^^j^Ys > 1) 
Ji 

+ [\\F'{7]s)\Ys + \F'{Ys)\x)^dxI{Ys < 1) 
Jys 
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<cE, 



(A.33) 



/■oo 

+ / {\Fm + C{l + Y,^^'))-dx 

Jl ^ 
f't / /"OO 1 rOO q 

j ds[AC{l + Y^+^) j -^dx + 2C j x^dx + -C 
+ 2C(1 - Ys)I{Ys < 1) + |F(0)| + C{1 + Y^+^)^ 



< cE„ 



ds 







^C + 5Cy;+2-^ + |F(0)| 
2 7 + i 



<c 



< oo 



since 7 + 2 G (0, 1) and where we also use the scahng property of X . (Note 
that without loss of generality we can assume that 7 is close enough to — 1 
so that 7 + 2 > 0.) Dividing both sides of (A. 32) by t, letting 1 1 and using 
the dominated convergence theorem, we get 



(A.34) 



[F{y - x) - F{y) + F'{y)xl{\x\ < l)]i^iidx) 



+ [FiO)-F{y)] i^iidx 



+ F'{y) / xl{\x\ < l)z^i(dx) 



for all y > 0. Splitting the integral over {—oo,y] into integrals over (—00, —y] 
and [—y, y] , noting that the third term of the resulting integral over (—00, —y] 
cancels with the final term in (A.34), it is easily seen using (A. 2) that the 
expression on the right-hand side of (A.34) coincides with the expression on 
the right-hand side of (A. 23). This completes the proof. □ 



REFERENCES 

[1] Abramowitz, M. and Stegun, I. A., eds. (1992). Handbook of Mathematical 
Functions with Formulas, Graphs, and Mathematical Tables. Dover, New York. 
MR1225604 

[2] Bernyk, v., Dalang, R. C. and Peskir, G. (2008). The law of the supremum 
of a stable Levy process with no negative jumps. Ann. Probab. 36 1777-1789. 
MR2440923 

[3] Bertoin, ,]. (1996). Levy Processes. Cambridge Tracts m Mathematics 121. Cam- 
bridge Univ. Press, Cambridge. MR1406564 



PREDICTING THE ULTIMATE SUPREMUM 



39 



[4] BOYCE, W. M. (1970). Stopping rules for selling bounds. Bell J. Econom. and Man- 
agement Sci. 1 27-53. MR0285259 

[5] DONEY, R. A. (2008). A note on the supremum of a stable process. Stochastics 80 
151-155. MR2402160 

[6] DU ToiT, J. and Peskir, G. (2007). The trap of complacency in predicting the 

maximum. Ann. Probab. 35 340-365. MR2303953 
[7] DU ToiT, J. and Peskir, G. (2008). Predicting the time of the ultimate maximum 
for Brownian motion with drift. In Mathematical Control Theory and Finance 
95-112. Springer, Berhn. MR2484106 
[8] DU ToiT, J., Peskir, G. and Shiryaev, A. N. (2008). Predicting the last zero of 

Brownian motion with drift. Stochastics 80 229-245. MR2402166 
[9] Gilbert, J. P. and Mosteller, F. (1966). Recognizing the maximum of a sequence. 
J. Amer. Statist. Assoc. 61 35-73. MR0198637 

[10] Gradshteyn, I. S. and Ryzhik, I. M. (1994). Table of Integrals, Series, and Prod- 
ucts. Academic Press, Boston, MA. MR1243179 

[11] Graversen, S. E., Peskir, G. and Shiryaev, A. N. (2001). Stopping Brownian mo- 
tion without anticipation as close as possible to its ultimate maximum. Theory 
Probab. Appl. 45 41-50. MR1810977 

[12] Griffeath, D. and Snell, J. L. (1974). Optimal stopping in the stock market. Ann. 
Probab. 2 1-13. MR0362766 

[13] Hochstadt, H. (1973). Integral Equations. Wiley, New York. MR0390680 

[14] Karlin, S. (1962). Stochastic models and optimal policy for selling an asset. In 
Studies in Applied Probability and Management Science 148-158. Stanford Univ. 
Press, Stanford, CA. MR0137594 

[15] Kyprianou, a. E. (2006). Introductory Lectures on Fluctuations of Levy Processes 
with Applications. Springer, Berlin. MR2250061 

[16] Patie, p. (2009). A few remarks on the supremum of stable processes. Statist. Probab. 
Lett. 79 1125-1128. MR2510779 

[17] Pedersen, J. L. (2003). Optimal prediction of the ultimate maximum of Brownian 
motion. Stoch. Stock. Rep. 75 205-219. MR1994906 

[18] Peskir, G. and Shiryaev, A. (2006). Optimal Stopping and Free-Boundary Prob- 
lems. Birkhauser, Basel. MR2256030 

[19] Protter, p. E. (2004). Stochastic Integration and Differential Equations. Springer, 
Berlin. 

[20] Revuz, D. and YOR, M. (1999). Continuous Martingales and Brownian Motion, 3rd 
ed. Grundlehren der Mathematischen Wissenschaften [Fundamental Principles 
of Mathematical Sciences] 293. Springer, Berlin. MR1725357 

[21] Sato, K.-i. (1999). Levy Processes and Infinitely Divisible Distributions. Cam- 
bridge Studies in Advanced Mathematics 68. Cambridge Univ. Press, Cambridge. 
MR1739520 

[22] Shiryaev, A. N. (2002). Quickest detection problems in the technical analysis of 

the financial data. In Mathematical Finance — Bachelier Congress, 2000 (Paris) 

487-521. Springer, Berlin. MR1960576 
[23] Shiryaev, A. N. (2009). On conditional-extremal problems of the quickest detection 

of nonpredictable times of the observable Brownian motion. Theory Probab. 

Appl. 53 663-678. 

[24] Urusov, M. a. (2005). On a property of the time of attaining the maximum by 
Brownian motion and some optimal stopping problems. Theory Probab. Appl. 
49 169-176. MR2141339 



40 



V. BERNYK, R. C. DALANG AND G. PESKIR 



[25] Watanabe, S. (1962). On stable processes with boundary conditions. J. Math. Soc. 
Japan 14 170-198. MR0144387 

[26] Williams, D. (1991). Probability with Martingales. Cambridge Univ. Press, Cam- 
bridge. MRl 155402 

V. Bernyk R. C. Dalang 

UBS AG Institut de Mathematiques 

EUROPASTRASSE 1 ECOLE POLYTECHNIQUE FeDERALE 

8152 Opfikon Station 8 

Switzerland 1015 Lausanne 

E-MAIL: violctta.bernyk@ubs.com Switzerland 

E-MAIL: robert.dalang@epfl.ch 

G. Peskir 

School of Mathematics 

The University of Manchester 

Oxford Road 

Manchester M13 9PL 

United Kingdom 

E-MAIL: goran@maths.man.ac.uk 



